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Abstract.  The  use  of  caged  dye  photo-activated  fluorophore  velocimetry  is  described  and 
representative  examples  are  presented.  After  a  brief  survey  of  recently  reported 
measurements,  a  more  detailed  example  of  the  flow  produced  in  a  cylinder  with  a  single 
rotating  end  wall  is  presented.  Simultaneous  stereoscopic  image  sets  in  the  (r,  z)  and  (r,  0) 
planes  have  been  obtained  over  a  Reynolds  number  range  of  roughly  102-105.  At  low 
Reynolds  numbers  (~0-2000),  the  steady,  axisymmetrical  flow  is  found  to  quantitatively 
agree  with  predictions  from  a  numerical  flow  solver.  At  higher  Reynolds  numbers  (from 
5  x  103  to  105),  the  flow  develops  considerable  turbulent  three-dimensional  structure. 
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1.  Introduction 

Recent  advances  in  laser  and  optical  technology  have 
motivated  the  development  of  new  quantitative  velocimetry 
techniques  based  on  two-  (and  three-)  dimensional  laser 
induced  fluorescence  (LIF)  imaging.  In  general,  these 
techniques  can  be  broadly  classified  as  either  ‘Doppler’  or 
‘time-of-flight’ .  Unlike  scattering  based  approaches,  such 
as  laser  Doppler  velocimetry,  Doppler  based  LIF  techniques 
most  commonly  utilize  the  spectral  shift  of  resonant  atomic 
or  molecular  absorption  features  that  are  associated  with 
the  component  of  fluid  motion  parallel  to  the  direction 
of  a  propagating  laser  sheet.  Since  such  features  are 
typically  somewhat  spectrally  broad  (of  order  1  GHz), 
reported  Doppler- shifted  LIF  measurements  have  mostly 
been  confined  to  high  speed  flows  (Zimmerman  and  Miles 
1980,  McDaniel  1983,  Paul  et  al  1989). 

Fluorescence  time-of-flight  techniques  determine  veloc¬ 
ity  by  imaging  the  displacement  elapsed  during  a  well  de¬ 
fined  time  interval.  Since,  in  principle,  the  time  interval  can 
be  made  arbitrarily  long,  time-of-flight  techniques  are  more 
suited  to  measuring  low  to  moderate  velocities,  such  as  those 
typical  of  incompressible  flows.  This  paper  will  focus  on 
what  has  been  termed  caged  dye  photo-activated  fluorophore 
(PAF)  flow  tagging  velocimetry.  Flow  tagging  refers  to  a 
class  of  time-of-flight  techniques  in  which  a  laser  beam  is 
used  to  ‘write’  a  pattern  into  a  flow  by  means  of  a  suitable 
optical  excitation.  After  an  appropriate  time  delay,  during 
which  time  the  original  pattern  is  advected  by  the  flow,  the 
new  position  is  ‘interrogated’  using  planar  laser  induced  fluo¬ 
rescence  (PLIF)  sheet  imaging.  Flow  tagging  based  on  vibra¬ 
tional  excitation  of  molecular  oxygen,  which  has  been  termed 
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Figure  1 .  A  schematic  illustration  of  the  use  of  caged  dye  PAFs 
for  flow  tagging  velocimetry. 


Raman  excitation  followed  by  laser  induced  electronic  fluo¬ 
rescence,  has  been  applied  both  to  low  speed  turbulent  air 
jets  (Noullez  et  al  1997)  and  to  supersonic  flow  (Miles  et  al 
1989).  Gas  phase  flow  tagging  based  on  in  situ  laser  produc¬ 
tion  of  hydroxyl  radicals  (Boedeker  1989,  Goss  et  al  1991) 
and,  more  recently,  ozone  (Ribarov  et  al  1999)  has  also  been 
reported.  An  inverse  type  of  tagging,  in  which  a  two-colour, 
multi-photon  absorption  is  used  to  photo-bleach  sodium  that 
is  seeded  into  a  flow,  has  also  been  reported  recently  (Barker 
et  al  1997). 

We  have  introduced  the  use  of  caged  dye  PAFs  as  flow 
tagging  tracers  in  water  and  alcohol  flows  (Lempert  et  al 
1995,  1998).  PAFs  are  nominally  fluorescent  dyes  that  have 
been  rendered  nonfluorescent  by  strategic  attachment  of  a 
chemical  ‘caging’  group.  The  chemical  caging  group  is 
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photolytically  cleaved  upon  absorption  of  a  single  photon 
of  ultraviolet  light.  After  photolysis,  the  fluorescent  dye 
is  recovered  and  can  be  tracked  indefinitely  using  ordinary 
laser  induced  fluorescence  imaging  techniques.  In  effect, 
as  illustrated  in  figure  1 ,  the  tagging  laser  is  used  to  photo- 
chemically  create  a  user  defined  pattern  of  fluorescent  dye. 
The  convection  of  the  ‘locally  seeded’  fluid  elements  is 
subsequently  interrogated  using  a  second  laser.  In  all  work 
reported  to  date,  the  tagging  has  been  performed  using  the 
third  harmonic  of  an  Nd:  YAG  laser  at  355  nm.  Interrogation 
has  utilized  visible  lasers,  most  commonly  argon-ion  lasers, 
flashlamp-pumped  dye  lasers,  or  the  second  harmonic  of 
Nd:YAG  (continuous  wave  or  pulsed). 

It  should  be  pointed  out  that  the  PAF  velocimetry 
technique  is  similar,  in  many  respects,  to  the  photochromic 
tracer  technique  (Popovich  and  Hummel  1967)  and  to  what 
has  been  referred  to  as  either  laser-induced  photochemical 
anemometry  (Falco  and  Nocera  1993)  or  molecular  tagging 
velocimetry  (Gendrich  etal  1997,  Koochesfahani  etal  1996). 

2.  Experimental  considerations 

The  majority  of  measurements  performed  to  date  using 
PAF  dyes  have  employed  some  form  of  caged  fluorescein. 
With  the  exception  of  the  addition  of  a  carboxy  dextran 
group  of  variable  relative  molecular  mass  and  one  (or  more) 
‘caging’  groups,  the  chemical  structure  of  the  uncaged  form  is 
identical  to  that  of  ordinary  fluorescein  dye.  Absorption  of  a 
single  UV  photon  results  in  permanent  photolytic  cleavage  of 
the  caging  groups.  Owing  to  its  high  brightness  and  fortuitous 
spectral  overlap  with  the  argon-ion  laser,  fluorescein  dye 
has  commonly  been  employed  for  passive  scalar  mixing 
and  flow  visualization  studies  (Koochesfahani  and  Dimotakis 
1985,  1986,  Dahm  and  Dimotakis  1990).  Like  fluorescein, 
the  uncaged  PAF  dye  has  a  fairly  broad  (~30  nm),  strong 
absorption  band  centred  at  480  nm  (blue)  and  exhibits  intense 
fluorescence  emission  (quantum  yield  ~0.90)  centred  at 
550  nm  (green).  More  recently  (Lempert  et  al  1998),  caged 
rhodamine  PAFs  with  optical  properties  similar  to  those  of 
caged  fluorescein,  but  with  their  optical  absorption  shifted  to 
~530  nm,  have  been  developed.  Caged  rhodamine  PAFs  are 
well  suited  to  interrogation  using  the  second  harmonic  of  a 
Nd:YAG  laser  at  532  nm. 

The  principal  advantage  of  caged  dye  PAF  tracers  is 
that  they  exhibit  very  strong  fluorescence.  As  discussed 
in  more  detail  by  Lempert  (2000),  the  radiative  decay  of 
the  upper  electronic  state  is  sufficiently  fast  (of  order  a 
few  nanoseconds)  that  many  photons  (102-108  or  more, 
depending  on  the  interrogation  source)  can  be  extracted  from 
each  uncaged  dye  molecule  during  a  single  interrogation. 
Additionally,  since  the  tagging  is  permanent,  PAFs  can 
be  used  for  measurements  in  exceedingly  slow  flows. 
For  example,  as  will  be  illustrated  in  section  3,  Harris 
et  al  (1996),  has  reported  quantitative  measurements  in 
electrohydrodynamic  flows  with  mean  velocities  of  order 
2-4  /im  s-1. 

There  are,  however,  some  significant  disadvantages 
associated  with  caged  dye  tracers.  Fundamentally,  the  most 
significant  is  the  finite  kinetic  rate  of  the  photochemical 
cage  breaking  step.  This  has  been  described  in  some  detail 


previously  (Lempert  et  al  1995,  1998)  and  has  the  effect 
of  imposing  a  minimum  delay  time  between  tagging  and 
interrogation.  To  date,  the  minimum  delay  time  which 
has  been  reported  in  caged  dye  PAF  studies  is  200  /zs 
(Lempert  et  al  1995).  This  corresponds  to  ~5%  of  the 
kinetic  1  /e  rise  time  of  the  fluorescence  signal.  The  second 
principal  disadvantage  of  caged  dye  PAFs  is  that  the  materials 
can  be  used  once  only,  since  the  uncaging  is  irreversible. 
This  is  compounded  by  the  high  cost  of  the  tracer  itself. 
PAFs  were  originally  developed  by  the  biological  sciences 
community  for  cell-related  studies  which  typically  require 
exceedingly  small  (of  order  a  milligram)  quantities  of  tracer. 
The  high  brightness  of  the  uncaged  dye  results  in  very  modest 
(~1  mg  l-1)  concentration  requirements,  but  the  cost  can 
still  be  prohibitive.  A  variety  of  caged  dye  PAFs  is  available 
commercially  from  Molecular  Probes,  Inc  (Eugene,  Oregon, 
USA). 

3.  A  survey  of  recent  measurements  employing 
caged  dye  PAFs 

Caged  dye  tracers  have  been  employed  in  a  wide  variety  of 
flow  problems,  with  characteristic  velocities  ranging  from  the 
order  of  micrometres  to  metres  per  second  and  characteristic 
length  scales  ranging  from  0.10  mm  to  1.0  m.  In  this  section 
we  give  a  very  brief  survey  of  measurements  which  have  been 
reported  to  date,  with  emphasis  on  providing  examples  of  the 
applicability  of  the  diagnostic  to  a  range  of  fluid  applications. 
More  detail  can  be  found  in  the  references  cited. 

At  the  low  end  of  the  length  scale,  Paul  et  al  (1998) 
have  obtained  quantitative  velocity  contours  in  pressure  and 
electrokinetically  driven  flow  through  open  capillaries  of 
order  100  /zm  diameter.  For  a  typical  pressure  driven  case, 
at  a  Reynolds  number  of  0.034,  a  classic  parabolic  profile 
was  obtained  with  a  centreline  velocity  of  450  /zm  s-1. 
Quantitative  velocity  data  were  obtained  by  recording  a  pair 
of  interrogation  images,  separated  by  20  ms,  and  using 
the  scalar  imaging  velocimetry  (SIV)  inversion  technique 
developed  by  Dahm  et  al  (1992).  The  spatial  resolution 
reported  was  approximately  5  /zm.  In  the  case  of 
electrokinetically  driven  flow,  it  was  noted  that  the  extraction 
of  quantitative  data  is  complicated  by  the  fact  that  caged 
fluorescein  is  charged  in  solution.  Separation  of  the  drift 
velocity  of  the  charged  tracer  from  the  true  fluid  convective 
velocity  was  found  to  be  a  considerable  problem.  The 
authors  also  noted  a  systematic  error  in  data  obtained  near  the 
capillary  wall,  which  also  was  attributed  to  possible  charge 
effects. 

Several  studies  of  internal  circulation  in  droplets  have 
been  reported.  Figure  2(a)  shows  a  tag/interrogate  image 
pair  obtained  from  a  free  falling  water  droplet  (Harris  et  al 
1996).  The  image  on  the  left-hand  side  of  figure  2(a)  is  due 
to  elastic  scattering  from  a  single  tagging  pulse.  The  image 
on  the  right-hand  side  of  figure  2(a)  shows  the  interrogated 
line  segment  29.5  ms  later.  The  droplet  in  figure  2  is 
approximately  5  mm  in  diameter  and  was  formed  from 
a  common  laboratory  burette  that  was  gravity  fed  from  a 
reservoir  containing  a  0.20  mg  l-1  (6.7  x  10~8  M)  solution 
of  caged  fluorescein  of  molecular  mass  3000.  As  described 
in  the  next  section,  quantitative  measurement  of  internal 
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Figure  2.  (a)  Tag  (left-hand  image)  and  interrogate  (right-hand 
image)  raw  image  pair  obtained  from  a  free  falling  water  droplet 
using  caged  fluorescein  PAF.  The  time  delay  is  29.5  ms.  (b)  The 
result  of  applying  ray  tracing  corrections  to  the  right-hand  image 
in  (a)  (from  Harris  et  al  (1996)). 

circulation  requires  that  the  raw  images  be  corrected  for 
the  influence  of  droplet  curvature.  In  effect,  the  droplet/air 
interface  acts  as  a  spherical  lens,  introducing  significant 
distortion.  Several  techniques  for  removing  this  optical 
distortion  have  been  reported,  including  rather  sophisticated 
approaches  which  remove  both  geometrical  and  intensity 
distortion  (Zhang  and  Melton  1994).  Application  of  a  simple 
ray-tracing  procedure  (Harris  et  al  1996),  summarized  in  the 
next  section,  gives  the  image  shown  in  figure  2(b).  It  can  be 
seen  that  the  region  very  close  to  the  droplet  edge  is  obscured 
by  total  internal  reflection.  Less  obviously,  the  magnitude 
(and  in  some  cases  the  sign)  of  the  measured  velocity  is 
greatly  effected.  Harris  et  al  (1996)  report  corrections  of  up 
to  1 00%.  Similar  data  have  been  obtained  for  alcohol  droplets 
(Lempert  et  al  1998)  using  a  caged  rhodamine  dye.  In  this 
case,  evaporation  of  alcohol  was  found  to  induce  exceedingly 
complex  internal  flow  patterns. 

As  an  example  of  an  exceedingly  low  flow  velocity,  a 
time  sequence  of  images  from  a  model  electrohydrodynamic 
flow  known  as  a  Taylor  cone  has  been  obtained  (Vizka  and 
Saville  1992).  The  cone  is  formed  by  immersing  a  pair  of 
electrodes  in  silicon  oil,  upon  the  bottom  of  which  is  placed 


a  water  droplet.  The  semispherical  droplet  is  ‘stretched’  into 
a  cone  by  application  of  an  approximately  2  kV  cm-1  field. 
The  fluid  velocity  reported  was  in  the  range  approximately 
2-4  /im  s-1. 

As  an  example  of  a  physically  larger  and  higher  Reynolds 
number  flow  field,  Biage  et  al  (1996)  obtained  a  large  data 
set  of  images  from  the  flow  produced  between  rotating 
concentric  cylinders  (Taylor-Couette  flow).  The  cylinder 
height  was  102  cm  and  the  inner  and  outer  cylinder  radii  were 
3. 14  and  8.18  cm,  respectively.  Sets  of  images  were  obtained 
over  a  Reynolds  number  range,  given  by  the  inner/outer 
cylinder  gap  spacing,  between  154  and  3.5  x  105  and  were 
used  to  obtain  instantaneous  velocity  profiles  and  spectral 
density  functions. 

As  a  second  example  of  measurement  in  turbulent  flow, 
Guilkey  et  al  (1996)  utilized  caged  dyes  to  study  scalar 
mixing  in  a  turbulent  pipe  at  a  Reynolds  number  of  7500. 
In  these  experiments,  a  cylindrical  optical  system  was  used 
to  create  a  spatial  ‘square- wave’  pattern  of  uncaged  dye.  The 
temporal  evolution  of  this  initial  pattern  was  subsequently 
monitored  at  several  downstream  locations.  It  was  found  that 
the  initial  square-wave  pattern  decayed  rapidly,  becoming 
dominated  by  higher  frequency  components  at  values  of  X/D 
equal  to  approximately  50.  At  X/D  =  71,  a  long  wavelength 
component  appeared,  as  predicted  in  previous  theoretical 
work  by  some  of  the  same  authors  (Kerstein  and  McMurty 
1994). 

As  a  final  example,  Dussaud  et  al  (1998)  utilized  caged 
dyes  to  visualize  the  spreading  of  involatile  and  volatile 
surface  films  on  water.  In  both  cases,  images  were  obtained 
by  writing  a  pair  of  vertical  lines  into  a  16  cm  diameter  tank 
filled  with  a  solution  of  caged  fluorescein.  A  single  droplet 
of  solvent  was  then  deposited  onto  the  quiescent  surface. 
Surface  tension  causes  the  droplet  to  rapidly  spread  into  a  thin 
surface  film.  The  sublayer  flow  pattern  induced  during  the 
spreading  was  captured  by  viewing  the  interrogated  image 
from  the  side.  For  the  involatile  fluid  (silicon  oil),  the 
experimental  velocity  data  agreed  well  with  Blasius  model 
predictions  for  a  laminar  boundary  layer.  For  a  volatile  fluid 
(p -xylene)  the  sublayer  flow  was  highly  rotational,  which  the 
authors  attributed  to  the  onset  of  a  thermal  instability  induced 
by  evaporation. 

4.  Rotating  end  wall  flow,  a  representative  example 

4.1.  Background 

In  this  section  we  use  the  representative  example  of  the  flow 
produced  in  a  cylinder  with  a  single  rotating  end  wall  to 
provide  a  more  detailed  look  at  caged  dye  flow  tagging.  This 
flow  is  attractive  because  it  provides  a  relatively  simple  model 
system  for  the  study  of  what  is  known  as  vortex  breakdown. 
Although  a  precise  definition  is  somewhat  elusive,  vortex 
breakdown  is  generally  characterized  by  the  formation  of 
a  stagnation  point,  followed  by  a  swelling  (or  ‘bursting’) 
of  the  vortex  and  the  creation  of  a  region  of  recirculation. 
After  Vogel’s  initial  observation  of  vortex  breakdown  in 
this  flow  (Vogel  1975),  Escudier  (1984)  utilized  fluorescent 
dye  visualization  to  systematically  explore  the  flow  over 
a  wide  range  of  the  parameter  space.  Recent  numerical 
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Figure  3.  A  simplified  schematic  diagram  of  the  lid  driven 
cylinder  flow  apparatus. 


simulations  by  Lopez  (1990)  and  analytical  predictions  by 
Brown  and  Lopez  (1990)  are  in  excellent  agreement  with  the 
flow  visualizations  of  Escudier  over  the  observed  Reynolds 
number  range  (<  3500).  The  purpose  of  the  current  studies 
is  to  assess  the  performance  of  the  caged  dye  technique  in 
such  a  well  characterized  model  environment  and  to  obtain 
new  quantitative  experimental  data,  particularly  at  higher 
Reynolds  numbers  for  which  the  flow  becomes  turbulent  and 
non-axisymmetrical. 


4.2.  Experimental  details 

Figure  3  is  a  simplified  schematic  diagram  of  the 
experimental  apparatus.  The  flow  is  contained  within  a 
precision  bore  quartz  cylinder  of  inner  diameter  6.6256  cm. 
The  cylinder  is  closed  on  the  bottom  with  a  stationary  quartz 
disc,  which  provides  optical  access  in  the  r-0  plane.  It  is 
closed  on  the  top  by  the  rotating  disc,  which  is  made  of 
Delrin  and  press  fitted  to  a  2.54  cm  precision  ground  stainless 
steel  shaft.  The  aspect  ratio,  8  =  H/R,  of  the  experiment 
is  adjusted  by  vertical  translation  of  the  drive  shaft,  which 
is  driven  with  a  dc  motor  via  a  belt  and  set  of  pulleys.  A 
combination  of  a  shaft  encoder,  computer  generated  reference 
frequency  and  closed  loop  feedback  control  ensures  that  the 
Reynolds  number  is  controlled  to  within  better  than  0.02%. 
The  working  fluid  for  all  data  presented  here  is  water,  with 
temperature  stabilized  by  a  water  jacket  and  measured  to 
an  accuracy  of  ±0.15  °C  by  precision  thermistors.  The 
combination  of  the  Reynolds  number  and  the  aspect  ratio 
completely  determines  the  flow  in  this  configuration. 

The  tagging  is  performed  by  focusing  a  single, 
approximately  20  mJ  pulse,  from  the  third-harmonic  output 
of  a  Nd:YAG  laser  operating  at  355  nm)  with  a  300  mm 


Figure  4.  An  interrogated  image  of  several  lines  tagged  at  1  s 
intervals  into  lid  driven  cylinder  flow  at  Re  =  960  and  8  =  2.0. 
Tagging  is  across  the  diameter  of  the  cylinder,  at  approximately 
one  quarter  of  the  cylinder  height,  measured  from  the  bottom. 

focal  length  lens.  The  interrogation  employs  a  flashlamp- 
pumped  dye  laser,  using  LD  490  laser  dye.  Typical  single 
pulse  energies  were  100  mJ  at  490  nm  with  a  pulse  duration 
of  1-2  fi s.  By  use  of  appropriate  optics,  the  interrogation 
beam  could  be  shaped  into  either  a  relatively  thin  sheet  or 
a  more  cylindrical  volume.  Sheet  optics  are  employed  at 
low  to  modest  Reynolds  number  ( Re  <  3000),  under  which 
conditions  the  flow  is  axisymmetrical.  At  higher  Reynolds 
number  (5  x  103  <  Re  <  105),  for  which  the  flow  becomes 
quite  three  dimensional,  the  cylindrical  interrogation  volume 
was  employed.  (This  is  the  case  illustrated  in  figure  3.)  A 
pair  of  ordinary  CCD  video  cameras  is  used  to  simultaneously 
capture  the  images  from  the  side  (the  r-z  plane)  and  bottom 
(the  r-0  plane).  A  simple  coloured  glass  filter  is  used  to 
attenuate  elastic  scattering  from  the  interrogation  laser,  while 
transmitting  more  than  80%  of  the  red-shifted  fluorescence. 
All  data  to  be  presented  were  obtained  using  dextran  DMNB- 
caged  fluorescein  dye  of  relative  molecular  mass  3000,  with  a 
concentration  of  approximately  1.6  mg  1_1 .  Complete  details 
of  the  experimental  apparatus  can  be  found  in  Harris  (1999). 

4.3.  Results  regarding  qualitative  visualization 

Figure  4  shows  a  series  of  lines  tagged  into  the  flow  at  1  s 
intervals  at  Re  =  960  and  8  =  2.0.  The  tagging  is  performed 
across  the  diameter  of  the  cylinder  at  the  position  of  the 
straight  horizontal  line  which  completely  spans  the  image 
at  a  point  located  approximately  one  quarter  of  the  way  from 
the  bottom  of  the  apparatus.  In  agreement  with  the  work  of 
Escudier,  no  stagnation  point  is  observed.  Rather,  the  tagged 
fluid  spirals  upwards,  rotating  out  of,  and  then  back  into, 
the  interrogation  sheet,  before  recirculating  back  downwards. 
Note  that,  after  approximately  30  s,  the  tagged  lines  remain 
well  defined,  indicating  a  lack  of  significant  mass  diffusion. 

Figure  5  shows  an  image  of  the  flow  at  a  Reynolds 
number  of  1890  and  <5  =  2.0.  In  this  case  a  strong  vortex 
breakdown  bubble  is  evident  near  the  bottom  and  a  fainter 
one  can  be  seen  at  the  approximate  vertical  midpoint.  The 
image  was  obtained  by  firing  approximately  60  tagging  shots 
into  the  flow  at  the  same  point  as  in  figure  4,  and  then  allowing 
the  flow  to  evolve  for  several  minutes.  Recirculation  keeps 
a  significant  fraction  of  the  uncaged  dye  constrained  to  the 
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Figure  5.  An  interrogated  image  of  lid  driven  cylinder  flow  at 
Re  =  1890  and  8  =  2.0.  Tagging  is  at  same  position  as  in  figure  4. 
Flow  was  tagged  several  times  and  then  allowed  to  evolve  for 
several  minutes  prior  to  interrogation. 


breakdown  bubble  zones.  As  the  Reynolds  number  is  further 
increased  (to  Re  —  3000),  the  flow  becomes  unsteady,  while 
maintaining  axisymmetry.  As  will  be  illustrated  in  the  next 
section,  at  even  higher  Reynolds  number  (Re  —  3500), 
the  flow  becomes  three  dimensional,  including  significant 
turbulent  fluctuations  in  the  boundary  layer. 


Figure  6.  An  interrogated  image  of  the  vertical  component  of 
displacement  in  lid  driven  cylinder  flow  at  Re  =  1410  and  <5  =  2.0. 
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Figure  7.  A  comparison  of  uncorrected  (+)  and  corrected,  ray 
traced  (0)  displacements  obtained  from  image  data  of  figure  6. 


4.4.  Quantitative  studies  on  steady  flow 

Extraction  of  quantitative  velocity  data  from  flow  tagging 
images  requires  the  consideration  of  several  issues,  the  most 
important  of  which  is  the  ambiguity  associated  with  motion 
parallel  to  the  direction  of  the  line.  This  has  been  discussed 
in  some  detail  by  several  authors  (Hill  and  Klewicki  1996, 
Koochesfahani  et  al  1996,  Paul  et  al  1998,  Harris  et  al 
1996)  and  has  resulted  in  several  strategies  which  will  only 
be  summarized  here.  The  simplest  approach  is  to  employ 
the  minimum  possible  time  delay,  At,  between  tagging  and 
interrogation,  since  the  uncertainty  in  the  displaced  line 
position,  E ,  is  given  by 

/  dw\  9 

E  =  (^)uAt  (1) 

where  w  is  the  component  of  velocity  perpendicular 
to  the  tagged  line  (typically  the  principal  flow  axis) 
and  u  is  the  velocity  component  parallel  to  the  line. 
Clearly,  minimizing  At  minimizes  this  contribution  to  the 
displacement  uncertainty.  The  finite  spatial  resolution  of 
cameras,  however,  causes  the  displacement  uncertainty  to 
increase  with  decreasing  At.  The  use  of  this  approach, 
therefore,  involves  an  inherent  trade-off  and  often  is  not 
a  satisfactory  solution.  A  second  approach  is  to  form  a 
more  extensive  grid  pattern,  utilizing  the  crossing  points  for 
unambiguous  two-component  velocity  determination.  This 
has  been  described  extensively  by  Kooschesfahani  and  co¬ 
workers  (Gendrich  and  Koochesfahani  1996,  Gendrich  et  al 
1997,  Koochesfahani  et  al  1996).  A  third  approach  is  to 
capture  a  sequence  of  images  and  employ  the  SIV  data 
inversion  procedure.  This,  as  mentioned  previously,  is 
described  in  detail  by  Dahm  et  al  (1992)  and  has  been 
employed  by  Paul  et  al  (1998). 


A  final  approach,  applicable  to  regimes  amenable  to 
direct  numerical  simulation  and  which  was  adopted  in  this 
portion  of  the  present  work,  is  direct  Lagrangian  frame 
comparison  of  the  observed  displacement  with  computational 
prediction.  This  is  achieved  by  ‘writing’  a  line  into  the 
computation  and  allowing  it  to  evolve  in  a  series  of  real¬ 
time  steps.  As  an  example,  figure  6  shows  a  representative 
image  obtained  from  the  side  viewing  camera  at  a  Reynolds 
number  of  1410  and  8  =  2.0.  The  line  was  tagged  across  the 
diameter  at  a  non-dimensionalized  vertical  position  equal  to 
1.15  R,  measured  from  the  rotating  lid.  The  time  between 
tagging  and  interrogation  was  0.40  s. 

The  displacement  at  each  lateral  pixel  location  along 
the  line  is  extracted  from  the  image  using  a  standard  least- 
squares  fitting  procedure  (Bevington  1969).  The  greyscale 
intensity  distribution  perpendicular  to  the  line  is  modelled 
as  a  Gaussian,  corresponding  to  the  idealized  tagging  laser 
intensity  distribution,  superimposed  on  a  linear  background. 
With  this  approach,  the  line  centre  can  typically  be  located 
with  a  minimum  statistical  uncertainty  of  ±0. 10  pixels.  More 
details,  including  examples  of  typical  raw  data  and  fits,  have 
been  presented  by  various  authors  (Lempert  et  al  1995, 
Lempert  2000,  Hill  and  Klewicki  1996,  Harris  1999). 

As  illustrated  in  the  previous  droplet  example,  when 
images  are  viewed  through  curved  surfaces  or  interfaces, 
the  effects  of  refraction  can  significantly  alter  the  observed 
displacement.  To  correct  for  this,  a  simple  pinhole  optical 
model  was  employed  (Harris  et  al  1996,  Harris  1999). 
The  basic  idea  is  to  computationally  ‘transfer’  the  intensity 
observed  at  the  CCD  image  plane  back  to  the  original  fluid 
object  plane.  Each  pixel  in  the  image  plane  is  mapped  onto  a 
ray  using  a  calibrated  lens  model.  This  ray  is  computationally 
propagated  through  a  model  of  the  experiment  until  an  optical 
surface  is  reached.  A  three-dimensional  refraction  angle 
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Figure  8.  A  Lagrangian  frame  comparison  between  ray  traced 
data  of  figure  7  (-)  and  a  numerical  computation  (O). 

is  then  calculated  using  Snell’s  law.  Linear  translation  is 
continued  until  the  ‘object’  plane  is  intersected.  If  only  a 
single  camera  is  available,  then  it  is  assumed  that  the  object 
plane  is  the  centre  of  the  interrogation  sheet.  This  implies  that 
the  line  was  tagged  in  the  centre  of  the  sheet  and  that  the  out- 
of-plane  motion  is  negligible.  If  two,  stereoscopic,  images 
are  available,  as  is  the  case  for  the  data  of  figure  6,  (a  bottom 
view,  not  shown,  was  also  obtained),  then  the  object  plane  is 
located  at  the  intersection  of  the  rays  corresponding  to  these 
two  views,  which  are  ray  traced  simultaneously.  Complete 
details  of  the  procedure,  including  bench  mark  tests,  can  be 
found  in  Harris  ( 1 999) .  Figure  7  shows  the  results  of  applying 
the  ray  tracing  procedure  to  the  data  from  figure  6.  Since  the 
curved  surface  is  a  cylinder,  the  principal  effect  is  to  distort 
the  radial  axis,  particularly  in  the  region  near  the  wall. 

Figure  8  shows  a  comparison  of  the  corrected  measured 
displacement  data  with  the  displacement  obtained  by 
simulating  the  motion  of  a  tagged  line  in  a  velocity 
field  obtained  by  solving  the  axisymmetrical  Navier-Stokes 
equations.  This  code  used  to  obtain  the  flow  field  is  very 
similar  to  one  described  by  Lopez  (1990).  In  figure  8, 
both  axes  have  been  non-dimensionalized  to  the  cylinder 
radius,  R ,  and  it  can  be  seen  that,  with  the  exception  of 
the  region  very  near  the  wall,  the  agreement  is  extremely 
good.  It  is  believed  that  the  discrepancy  near  the  wall  is 
entirely  due  to  errors  associated  with  ray  tracing.  More 
detailed  examination  (Harris  1999)  shows  that  the  average 
deviation  between  measured  and  computed  displacement  is 
approximately  0.00147?. 

4.5.  Results  for  the  unsteady  regime 


Figure  9.  Seven  consecutive  images  of  azimuthal  motion  (in  the 
r-0  plane)  of  a  line  tagged  into  lid  driven  cylinder  flow  at  four 
different  Reynolds  numbers.  Re  =  5  x  103  (a),  104  (b), 

5  x  104  (c)  and  105  (d). 


magnitude  of  large  scale  unsteadiness  has  increased  and 
small  scale  features  are  beginning  to  evolve  near  the  walls. 
At  Re  =  5  x  104  (figure  9(c)),  small  scale  features  are 
quite  evident.  At  Re  =  105  (figure  9(d)),  the  scales  of 
the  fine  structure  are  approaching  the  width  of  the  tagged 
fluid  element.  It  should  be  noted  that  the  increased  noise  in 
figure  9(d)  is  a  result  of  scattering  from  small  air  bubbles 
aspirated  into  the  flow  from  above  the  rotating  end  wall. 
Measurements  at  the  highest  Reynolds  numbers  employed  a 
time  delay  between  tagging  and  interrogation  of  less  than  the 
1/60  s  video  framing  interval.  In  this  case,  scattering  from  the 
small  bubbles  by  the  tagging  laser  pulse  is  superimposed  on 
the  interrogation  image,  despite  the  use  of  the  coloured  glass 
filter.  These  bubbles  were  produced  only  when  the  apparatus 
was  run  at  its  highest  Reynolds  number  setting  and  do  not 
appear  to  be  significant  enough  to  perturb  the  flow. 

The  two-point  turbulent  displacement  correlation 
function,  C(/y,  rj),  provides  a  convenient  way  to  present 
quantitative  information  about  the  structure  of  this  flow.  The 
correlation  function  is  defined  as 


Figure  9  shows  seven  consecutive  images  of  azimuthal 
motion  (in  the  r-0  plane)  obtained  from  the  bottom  camera  at 
four  different  Reynolds  numbers  ranging  from  5  x  103  to  105. 
Although  it  is  not  readily  apparent  from  figure  9,  even  at  the 
lowest  Reynolds  number  of  5  x  103  the  flow  has  developed 
significant  unsteadiness.  (This  is  much  more  apparent  when 
image  sequences  are  displayed  in  rapid  succession  on  a 
screen.)  There  are,  however,  no  discernible  small  scale 
features.  At  a  Reynolds  number  of  104  (figure  9(b))  the 


x(rj)x(rj) 

[x(ri)*x(rj)*]W 


(2) 


where  x(rt)  is  the  displacement  at  the  ith  radial  coordinate 
and  the  overbars  indicate  an  ensemble  average  over  the 
individual  images  obtained  for  each  flow  condition  (150 
frames  at  Reynolds  number  104,  ~200  for  Reynolds  numbers 
5  x  104  and  105).  As  the  time  delay  between  the  tagging  and 
interrogation  steps  approaches  zero,  the  value  of  C(/y,  rj) 
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Figure  10.  Spatial  correlation  plots  obtained  from  ensembles  of 
instantaneous  images  similar  to  those  of  figure  9.  At  Re  =  104 
(top),  5  x  104  (middle)  and  105  (bottom).  Red  =  correlation  of 
1.0,  blue  =  correlation  of  —  1 .0. 


approaches  that  of  the  familiar  two-point  turbulent  velocity 
correlation.  Figure  10  shows  colour  plots  of  the  correlation 
function  of  the  azimuthal  displacement  as  a  function  of 
the  non-dimensionalized  radial  coordinate  at  three  Reynolds 
numbers.  The  red  diagonal  bands  are  a  result  of  the 
strongly  correlated  motion  of  nearby  points.  The  perfect 
reflection  symmetry  across  the  diagonal  band  is  due  to  the 
property  of  the  correlation  function  C(r/,ry)  =  C(ry-,  rt). 
The  weaker  symmetry,  C (rt ,  rj )  =  C(— rt,  —  rj)  is  an 
indication  of  the  axisymmetry  of  the  flow  itself.  The 
top  plot  of  figure  10  corresponds  to  a  Reynolds  number 
of  104.  Of  particular  interest  are  the  question-mark- 
shaped  features  in  the  upper  right-hand  and  bottom  left- 
hand  corners,  indicating  that  the  turbulent  fluctuations  on 
opposite  sides  of  the  cylinder  are  correlated  and  in  phase.  At 
Re  =  5  x  104  (middle  image),  the  width  of  the  strongly 


correlated  diagonal  is  decreased,  indicating  a  decrease  in 
the  characteristic  turbulent  length  scale.  There  are  also 
multiple  off-diagonal  features,  indicating  that  the  flow  is 
developing  some  rather  complicated  large  scale  structure. 
The  bottom  plot  of  figure  10  corresponds  to  a  Reynolds 
number  of  105  and  exhibits  regions  of  significant  positive 
and  negative  correlation.  Negative  values  correspond  to 
turbulent  fluctuations  which  are  correlated,  but  are  180°  out 
of  phase.  For  example,  negative  correlation  at  points  such 
as  (ji  =  —0.25,  rj  =  0.25)  indicate  that,  at  this  Reynolds 
number,  the  motions  of  structures  on  opposite  sides  of  the 
core  are  coupled  together  in  such  a  way  as  to  produce  out-of- 
phase  displacement. 

5.  Conclusions 

Caged  dye  flow  tagging  velocimetry  has  been  applied  to 
a  wide  variety  of  flow  problems,  spanning  a  large  range 
in  characteristic  length  and  time  scales.  Quantitative 
velocity  data  for  microchannel  flow,  free  falling  and 
electrohydrodynamic  droplets  and  rotating  flows  spanning 
a  Reynolds  number  range  from  of  order  0.10  to  of  order 
105  have  been  reported.  Visualizations  of  the  highly  vortical 
sublayer  flow  produced  by  an  evaporating  droplet  spreading 
over  a  surface  have  also  been  reported,  as  has  a  novel 
approach  for  measuring  scalar  mixing  in  turbulent  pipe 
flow. 

As  a  more  detailed  example,  caged  dye  flow  tagging 
has  been  used  to  study  the  complex  flow  produced  in  a 
cylinder  with  a  single  rotating  end  wall.  At  low  Reynolds 
number  (of  order  100-2000),  visualizations  in  the  r-z 
plane  are  in  agreement  with  the  dye  studies  of  Escudier. 
Simultaneous  stereoscopic  images  in  the  r-z  and  r-6  planes 
have  been  obtained  under  these  conditions.  After  removal  of 
geometrical  distortion  produced  by  the  curved  walls  of  the 
container,  quantitative  agreement  between  the  experimental 
images  and  those  predicted  using  a  numerical  flow  solver 
is  found.  At  still  higher  Reynolds  number  (5  x  103  to  105), 
the  flow  develops  considerable  smaller  scale,  turbulent  three- 
dimensional  structure,  particularly  in  the  boundary  layer  on 
the  cylinder  walls.  Sets  of  images  have  been  obtained  and 
used  to  determine  statistical  quantities  such  as  the  spatial 
correlation  function. 
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